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Abstract 

We investigate diverse relations of the colour-flavour transformations (CFT) introduced 
by Zirnbauer in [41, 42] to various topics in random matrix theory and multivariate analysis, 
such as measures on truncations of unitary random matrices, Jacobi ensembles of random 
matrices, Berezin reproducing kernels and a generalization of the Selberg integral due to 
Kaneko, Kadell and Yan involving the Schur functions. Apart from suggesting explicit 
formulas for bosonic CFT for the unitary group in the range of parameters beyond that in 
[42] , we also suggest an alternative variant of the transformation, with integration going over 
an unbounded domain of a pair of Hcrmitian matrices. The latter makes possible evaluation 
of certain averages in random matrix theory. 
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1 Introduction 

Let (xj)^!-! and (y^JLi be two sets of column- vectors in C N . By convention, the column- 
vectors are regarded as matrices consisting of a single column, and we shall use the star to 
denote the Hermitian conjugate (complex conjugate transpose) of a matrix. Let U (N) stand for 
the group of unitary matrices of size N x N equipped with the Haar measure dfiniU) fixed by 
the normalization condition Ju^ N ^ dfj,n(U) = 1. Further, consider the matrix ball Q*Q < I m in 
C mxm , the space of complex m x m matrices, equipped with the unit mass measure 

d/4, m (Q) = const. det(J m - Q*Q) N - 2m (dQ), N > 2m, (1) 
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where I m is identity matrix and (dQ) is the cartesian volume element in C mxm , 

m 

(dQ) = ] [ dReQj k dlmQj k . 

j,k=i 

With these notations in mind, the following remarkable identity, called the bosonic Colour- 
Flavour Transformation (bCFT), 

f e Z?=M u *i+*S u 'Vi) dnHiP) = [ e £&=i«*^+«-)*ito) d^ m {Q), (2) 

Ju(N) JQ*Q<I m 

is known to hold for N > 2m. Similarly, the fermionic Colour-Flavour transformation (fCFT) 
asserts that 

f e Er=i(*W+^*x,) dfiH{u) = f e E£ fc=1 (Q*xto-«-)*^-) d M £ m (Q), (3) 

JU(N) Jc mxm 

where now Xji ^ji Xj an d rfj are vectors with anti-commuting components. The Xj an d ipj 
are column- vectors and the Xj an d iftj are row- vectors. In contrast to (2), there is no restriction 
on m and N in (3). The integration on the right-hand side in (3) is over the entire space of 
complex mxm matrices and 

d^N,m(Q) = const. det- N - 2m (I m + Q*Q) (dQ). (4) 

Both identities, together with the unifying supersymmetric variant of Colour-Flavour Trans- 
formation (CFT) and extensions to other classical groups, were originally discovered by Zirn- 
bauer in 1996 [41, 42] and proved by a skillful use of the machinery of generalized coherent 
states [33]. Following Zirnbauer's approach, variants of the CFT were also obtained for the 
special unitary group in [7, 34, 39]. Later on it was realized, again by Zirnbauer, see e.g [44, 10], 
that actually all forms of the CFT are just manifestations of a very deep algebraic fact related 
to the so-called Howe duality [20]. Since their introduction, the Colour-Flavour Transformations 
proved to be a very useful tool, finding diverse applications in such areas of physics as lattice 
gauge theory [42, 7, 36], random network models [1, 43], quantum chaos models [2, 3], and the 
random matrix theory [41, 16]. 

In this paper we revisit the simplest case of the unitary group. Our inspiration comes from 
noticing a certain similarity between a few results known in the random matrix theory and the 
bosonic and fermionic versions of the CFT. The central role in our investigation is played by 
certain integrals involving the so-called Schur functions s\(A). The latter are explicitly defined 
for any mxm matrix A in terms of its eigenvalues a±, . . . ,a m as 

, , / m+Xj—j\m 

det [a { jj 7=1 

s x (A) = s x (a 1 ,...,a m ) = — — - • m ' 3 , Xi > A 2 > . . . > A m > 0, (5) 

det (a,- ) ■ ■ , 

with A being a partition, i.e. a non-increasing sequence of non-negative integers Xj. The Schur 
functions s\(A) are symmetric polynomials in the eigenvalues of A = (Aij), and are also poly- 
nomials in A^. A concise introduction to the theory of symmetric functions can be found in the 
first chapter of [27], see also [8]. 
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The usefulness of Schur functions for our purposes can be traced to the fact that they are 
characters of irreducible representations of the general linear group and its unitary subgroup 
and, as a consequence, possess certain properties of orthogonality, see, in particular, equations 
(51) and (52). This makes Schur function expansions a powerful tool for evaluating integrals 
over unitary groups as has been already demonstrated in [26, 4, 36, 39, 31, 32, 19]. In this paper, 
we use the Schur function expansion technique to extend the bosonic version of CFT (2) to the 
range N < 2m < 2N. We also reveal the relation of bCFT and fCFT to several important matrix 
integrals due to Berezin, and also to certain generalizations of the famous Selberg integral due 
to Kaneko, Kadell and Yan. Finally, we also derive a new variant of the bosonic CFT, replacing 
the integration on the right-hand side in (2) with one going over an unbounded matrix domain 
parameterized by a pair of Hermitian matrices. Such a representation is useful for studying 
regularized inverse spectral determinants of complex random matrices, the subject of our earlier 
work [16]. The latter work has in fact quite a few points of intersection with some of the topics 
discussed in the present paper. 

The fact that the Schur functions play central role in our way of understanding and extending 
both bCFT and fCFT can be traced back to the above mentioned Howe duality, although we do 
not exploit the latter explicitly in the present paper. Without going into any detailed discussion, 
we would like to mention that in one of its incarnations Howe duality can be looked at, see e.g. 
[9] or Chapter 43 in book [8], as the ultimate reason for the validity of the so-called Cauchy 
identities 



m n ^ 

mi 

i=i j=i j 



^s A (ti,...,t m )s A (xi,...,x n ) (6) 



[]n( 1 + f!;i; i) = ^2 s x(h,---,t m )s X '{x 1 ,...,x n ), (7) 

i=l j=l x 

where A' stands for the conjugate partition. Being "a point of very direct connection between 
representation theory and combinatorics" [8], it is therefore no surprise that both identities (6) 
and (7) play an important role in obtaining the results of this paper. 

The organization of the paper is as follows. In the next section we will provide, for the 
reader's convenience, an overview and short discussion of the most important results following 
from our approach to CFT. The subsequent sections will be devoted to technical derivations of 
the main formulas. 



2 An overview of the main results. 

Our first observation relates the bCFT to truncations of random unitary matrices. Truncations 
of random unitary matrices emerged recently in a different context of quantum chaotic scattering. 
Various matrix distributions arising from such truncations were the subject of a few recent works, 
see [45, 13, 37, 11]. 

The relation between the bCFT and truncations of unitary matrices becomes more apparent 
if one writes (2) using matrix notation in the exponential. Introducing two N x m matrices X 
and Y with columns x±, . . . , x m and y 1 , . . . , y m respectively, one can write 

m 

^(ypxj + x*U*yj) = Ti{Y*UX + X*U*Y) = Tr \XY*U + U*YX*) 

3=1 
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and 

X*X = (x* j x k )% k=1 and Y*Y = (y*y k )™ k=1 . 
Then the bCFT (2) takes the following form 



L 



e Tr(XY*U+U*YX*) ^ ([/) = j ^ XQ + Q*Y*Y) ^^(Q). (8) 



U(N) JQ*Q<I, 

Note that the left-hand side in (8) is well defined for any m whilst the integration measure on 
the right-hand side is singular if 2m > N and care must be taken if one wants to interpret 
the above formula for 2m > N. We argue in Section 3 that for m < N one can rewrite the 
right-hand side in a form free from singularities. Indeed, the matrix XY* and its Hermitian 
conjugate Y*X have rank m < N, or, to be more precise, at most m. In view of the invariance of 
the Haar measure, this means that the integral on the left-hand side in (8) goes effectively over 
the principal (top left) mxm sub-block Q of the unitary matrix U. With this observation in 
hand, a straightforward application of the Schur function expansion for the exponential function 
exp Tr M allows one to recast the bCFT in the following form 

[ e Tr ^+^) dp H (U) = [ e Tr(X*XQ +Q *Y*Y) dp N , mxm {Q), (9) 

JU(N) JQ*Q<Im 

where dpM,mxm{Q) is the image of the Haar measure under the truncation U i— > Q. It is known 
from [13] and [11], see also [30], that dp N)inxm (Q) oc det(7 m - Q*Q) N ~ 2m (dQ) for N > 2m. 
Thus, as one would expect, (9) reverts back to the original version of the bCFT in the interval 
N > 2m. 

Identity (9) which holds for any m < N (note that dpN, m -xm{Q) is a unit mass measure and 
is free from singularities) provides the basis for our extension of the bCFT to the range m < N . 
The boundary case m = N is straightforward. Indeed, if m = N then no truncation is involved 
and (9) takes the form 

j e Tr(Y*UX + X*U*Y) dtiH{u) = f ^X'XU+U'Y'Y)^^ 
Ju(N) Ju(N) 

The intermediate range TV" < 2m < 2N requires additional calculations which are of interest 
on their own. Though in this range the measure dpN,mxm(Q) can still be found explicitly, the 
resulting expression is too complicated to be practically used. However, for our purposes we 
only need to know the radial part of dpN,mxm{Q) and this latter can be found in the explicit 
form and then used to rewrite the right-hand side of (9). 

As we discuss in Section 3, in the range iV < 2m < 2N the measure dpN, m xm(Q) is supported 
by the set T on the boundary of the matrix ball Q*Q < I m defined by the condition that 
rank(/ m — Q*Q) = N — m. This set can be parametrized by the matrices 

where Z runs through the matrix ball Z*Z < In-th in C^ N ~ m ^ x( - N ~ m \ the space of complex 
(N — m) x (N — m) matrices, and U and V run through the unitary group U (m). However, this 
parametrization is not one-to-one: different Quzv* m &y represent the same point Q in T. In 
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the terminology of [21] (page 118) the set of matrices Quzv* is a covering of T. Exploiting such 
a parametrization, we obtain the following variant of the bCFT in the range ./V < 2m < 2N 

I e T ^ Y * ux+x * u ^dp H (U) = (11) 

Ju(N) 

III e^^v^uzv^y) d^ N _ m (Z)dMUW H (V), 

JU(m) JU(m) JZ*Z<I N _ m 

where Quzv* is as defined in (10) and dji^ N _ m (Z) is the unit mass measure on the matrix ball 
Z*Z < In-wl defined in (1), 

d^N, N -m( z ) = con st- det(/7V-m - Z*Z) 2m - N (dZ), N <2m< 2N. (12) 
It has to be mentioned that another variant of extension of the bCFT to the range m < N, 

I e T « Y * ux+x * u * Y Up H (U) = [ det(VY*X) m - N e^ x * XVf+VY ^d m (V), (13) 

Ju(N) Ju(m) 

was obtained recently in [39]. Clearly, apart from the boundary case of N = m this formula is 
different from ours. Interestingly, as was observed in [39], in the range m < N, the integral on 
the left-hand side in (13) does not change if the integration over the unitary group is replaced 
by the integration over the special unitary group: 

f e T < Y * ux+x * u * Y Up H (U)= f e^ Y * ux+x * u * Y U m (U), m < N. (14) 

Ju(N) JSU(N) 

This means that for m < N our formula (11) holds without further changes if one replaces U(N) 
by SU(N) in the integral on the left-hand side. 

The original CFT (2) and its extension (9) almost hide the fact that both integrals, the one 
on the left-hand side and the one on the right hand side, depend only on the eigenvalues of 
the matrix product X*XY*Y. In fact, the bosonic CFT can be stated in the following, slightly 
more abstract form. If A and B are two N x N matrices of rank m < N then 
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e Tr(AU + U*B) dMu) = f e Tr(CQ+Q*D) dp N , mxm {Q) (15) 

JQ*Q<I m 



U{N) JQ*Q<1, 

for any pair of m x m matrices C and D such that the eigenvalues of CD coincide with the 
non-zero eigenvalues of AB. In the range N < 2m < 2N the integration over matrices Q on the 
right-hand side can be replaced by integration over the matrices Quzv* as in (11). 

Our next observation is that the CFT is related to several interesting matrix integrals. 

It will be apparent from our calculations in Section 3 that the bosonic CFT (9) implies the 
following identity for the Schur functions 

I s x (Q*Q) d PN , mxm (Q) = m<N, (16) 

JQ*Q<I m s x {1n) 

and vice versa, (16) implies (9). In fact, (16) is a particular case (corresponding to m = n) of a 
more general relation 

/ s x (Q*Q) dp N , nxm (Q) = s ^ Im >^ In \ m ,n<N. (17) 

jQ*Q<I m sx(1n) 
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Here the matrices Q are n x m and dpN,nxm(Q) is the image of the Haar measure under the 
truncation of unitary matrix to its principal n x m sub-block, U i— ► Q. This relation and its 
generalization 

f s x {XQYQ*) dp N , nxm (Q) = S ^ X >^\ m,n<N, (18) 

jQ*Q<Im S\(1 N ) 

are simple corollaries of the invariance of the measure dpN, n y.m{Q) with respect to the right 
and left multiplication by unitary matrices. Another corollary of this invariance (and of the 
orthogonality of the Schur functions, see (52)) are the following orthogonality relations 



Ic 



s\(M*L) 



s\(LQ)s fjL {MQ) dp N 

,nxm (Q)=Sx,v \ T m,n<N, (19) 

Q*Q<I m ' Sx(lN) 

which hold for arbitrary mxn matrices L and M. Identities (18) and (19) are derived in Section 
3. 

We show that the integration formulas (17)— (19) imply several non-trivial matrix integrals, 
some of which we believe to be new. In particular, by making use of the Schur function expansion 
(6), we derive the following identity 



L 



d»H(U) _ 



U(N) det(I N - AU) m det(/jv - U*B*) n 

— L<_^ w-vzob-a) det(/ " z " z)N -"- n {dZ) 

which reduces the group integral on the left to an average over the Jacobi ensemble of random 
matrices Z of size k x k, k = min(m,n). Identity (20) holds for N > n + m and A* A < In, 
B*B < In and generalizes our earlier result [16] from n = m to the case n/m. A similar 
identity holds in the range N < m + n < 2N. 

In [16] we obtained an identity which is dual to (16), 

/ s x (Z*Z) d^ m (Z) = %1. (21) 

Here dp N m , see (4) , is the measure which appears on the right-hand side in the fermionic version 
of the CFT (3). By making use of the Schur function expansion (7), identity (21) is equivalent 
to the matrix integral 

f det(7jv + AU) m det(/jv + U*B*) m dp H (U) = [ det(7 + Z*Z B*A) dp Nm (Z) (22) 

Ju(N) JC mxm 

which is dual to the m = n version of (20). The emergence of dp Nm (Z) on the right-hand side 
in (22) is not coincidental. In fact, we show in Section 4 that the fermionic version of the CFT 
(3) implies (22) directly. We expect that identity (21) (or, equivalently, the matrix integral (22)) 
should in turn imply the fermionic CFT, by analogy to the relation between the bosonic CFT 
and identity (16). Unfortunately, we have succeeded in verifying such equivalence only for the 
simplest case m = 1. 
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In the same way as identity (16) allows for an extension to rectangular matrices, identity 
(21) allows for a similar extension: 

/ sx(Q*Q) d^ nxm (Q) = S ^>f™\ (23) 
where dpfj nX m(Q) * s um * m &ss measure 

df*N,nxm(Q) = const - det(I m + Q*Q)- N - m - n (dQ) 

on the space of complex n x m matrices. Identity (23) holds for any positive integers N, m and 
n, and in turn, implies the identity dual to (20): 

/ det(Jjv + AU) m det(/jv + U*B*) n dp H (U) = (24) 

Ju(N) 

f r\pf( 7* 7\\ n ~ m \ 

const - / det ( J + z * z ® B * A ) a , ; 7 jw +m+w 

J C min(m,»)xmin(m, n ) det(i + Z*Z)" +m+n 

again reducing evaluation of the integral over the unitary group on the left to evaluation of an 
integral over a Jacobi ensemble of random matrices. This is a generalization of our earlier result 
from [16]. 

Identities (16) and (21) yield another pair of matrix integrals, again by the way of the Schur 
function expansions (6)-(7): 



So 



dpN,nxm(Q) _ I 



l Q * Q < Im det(/„ - Z X Q*Y det(J n - QZ*)" det(J n - Z X Z\Y ' 
where Z*Zi < I m , and 

det(/ m + Z X Q*) N det(/ m + QZ*) N d^ nxm (Q) = det(7 m + Z X Z*) N . (26) 



/ 

Jo 



These matrix integrals are not new. They are a variant of integrals obtained by Berezin in his 
work on quantization in complex symmetric spaces [5, 6]. Berezin proved (26) for integer N and 
(26) for a range of real N that includes N >n + m. In Section 4 we discuss this link and quote 
some of Berezin's results. 

Since identities (17) and (23) are so useful in the context of Schur function expansions, we 
think it is worth to have a closer look at them. Without loss of generality we can assume 
that n > m. In the bosonic case we also assume that N > n + m, so that the integration 
measure dpN^ nxm {Q) in (17) is replaced by const. det(7 m — Q*Q) N ~ m ~ n (dQ). In the range 
N < n + m < 2N one can obtain slightly different formulas by using parametrization (10), see 
especially the integration formula (50). 

The integration in (17) and (23) is effectively over the eigenvalues of Q*Q. By making the 
corresponding change of variables (see, [21] or [12]) one brings the matrix integral in (17) to 



-k-f 1 - C »a(*i, . . . , x m ) n *r m (i - x ^ N ' m ~ n n - n ^ = sx{l :): x{ i rn) 

C n,mJ0 JO j=1 l<i<j<m j=l S X (l N ) 

(27) 
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and the one in (23) to 



j roc roo «<■ X S A (1 )s\(l ) 

— ... SX ( X1 ,..., Xm )ii 3 N+m+n n (^- x j) 2 n^= \ s i(i^ m ■ 

V J ° - 70 j=l y l+X 3> l<i<j<m j=l Sy{l N ) 

(28) 

The normalization constants in these formulas can be computed with the help of the celebrated 
Selberg integral, and are given in (87) and (88). 

Interestingly, the integral in (27) is the 7 = 1 case of the following extension of the Selberg 
integral due to Kaneko [25] (see also related works by Kadell [24] and Yan [40]) 



/•••/ J~x{xi, ■ ■ ■ ,Xm) a^ -1 (l - Xj)^ 1 II \xi - Xj\ 21 Y[ dx J 

^° ^° .7=1 l<i<j<m .j=l 

x m 



(29) 



j=i i<*<i<™ 3=1 

F(ij + l)r(Aj + p + 7(771 — «))r(<7 + 7(771 — i)) 



r(l + 7 )r(Ai + p + g + 7(2^ - i - 1)) 

Here (2) are the Jack symmetric functions [38], (x) are proportional to the Schur functions 
if 7 = 1, Jl(x) = H\s\(x), where H\ is a coefficient independent of x. Similarly, (28) is the 
7 = 1 case of the integral 



/•oo /-oo 1 m jP- 1 m 

... jz( tl ,...,t m )ii ; +g+2(m _ 1)7 n i^-^i 27 

J0 J{J .7=1 v T J 7 l<«<i<m 3=1 



To the best of our knowledge, the latter integral is not evaluated yet for 7 ^ 1 and we consider 
its computation as an interesting open problem. In particular, we note that although for the 
zero partition Ai = A2 = . • • = A m = the substitution x = 1/(1 — t) reduces the integral in 
(30) to the one in (29), this substitution does not preserve the Jack symmetric functions, and 
hence the integral in (30) for non-zero A requires a separate evaluation. 

Surprisingly, the 7 = 1 case is simple in the sense that the integral in (30) (and the one in 
(29)) can be evaluated by elementary means, as we demonstrate in the end of Section 4. 

Section 5 of our paper is devoted to yet another alternative ("deformed" ) version of the 
bosonic CFT with a different integration manifold in the integral on the right-hand side of (2). 
Consider the manifold of matrices (Qi,^) in C mxm x C mxm parametrized as follows 

Qi = TPT*, Q 2 = (T*)- 1 PT -1 , (31) 

where T runs through the general linear group, T 6 GL m (C), and P runs through the set of 
diagonal matrices 

P = diag(pi, . . . ,p m ); -1 < Pj < 1, j = 1, . . . ,m, (32) 
and introduce the integration "measure" 

m 

(dQ 1 dQ 2 )= H (pZ-plfUpjdpjdfiHiT), (33) 

l<j<k<m 3=1 

where d^niT) = (dT)/ det(T*T) m is the invariant measure on GL m (C). Note that pj change 
sign, and therefore (dQidQ?) is not a proper positive measure, see [15] for a discussion. Then 
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we claim the validity of the following variant of the bosonic CFT (N > 2m) 

I e-^ Y * ux+x * u * Y Up H (U) = [ det(/ m - Q 1 Q 2 ) N - 2m e- i ^ Y ' Y+ ^ x * x \dQ 1 dQ 2 ). 

Ju(N) J 

(34) 

Note the most important features of this formula are: (i) integration domain in (34) is un- 
bounded, in contrast to the bounded domain Q*Q < I m as in the standard bCFT (8) and (ii) 
the matrices Q\ and Q 2 are Hermitian. We finish that section by demonstrating that the last 
property makes the new version of bCFT indispensable when evaluating expectations values of 
negative powers of certain (regularized) spectral determinants. 



3 bCFT from truncations of random unitary matrices 

In this section we first address the problem of evaluation of integrals of the form 

f(Q*Q) dp N>nxm (Q), (35) 



Q*Q<I m 

where dpM,nxm{Q) is the image of the Haar measure under the truncation of unitary matrices 
as defined below. Then we apply the obtained formulas to derive the bCFT in the form of 
equations (9) and (11). 

Let U be an N x N unitary matrix and m and n be positive integer numbers, m < n < N . 
Partition U into the four blocks 

u -(ps) < 36 > 

where the top left block Q is n x m. Partition (36) defines the map oj : U — > Q from the unitary 
group into the matrix ball Q*Q < I m in C nxm , the space of complex n x m matrices. Under 
this map the Haar measure dpn(U) on U(N) induces a measure on the matrix ball Q*Q < I m 
which we shall denote by dpN,nxm(Q)- The unitarity of U imposes constraints on its sub-blocks. 
In particular, 

Q*Q + P *p = Im . (37) 

If N > m + n then, generically, the matrix P*P has rank m, and the image of U(N) under 
the map u is the entire matrix ball Q*Q < I m . In this case the measure dpM,nxm{Q) has been 
previously computed in [30, 13] for square matrices Q and in [11] for rectangular matrices, 

dpN,nxm(Q) = const. det(/ m - Q*Q) N ~ m - n (dQ), N>n + m, (38) 

where (dQ) is the cartesian volume element in C nxm . By making an appropriate change of 
variables, see e.g. [12] or one of the calculations below, the integral in (35) can then be reduced 
to the familiar form 



f f(Z*Z) det(Z*Z) n ~ rn det(7 m - Z* Z) N ~" {dZ) (39) 

Jz*z<i m 



of an average over the so-called Jacobi ensemble [12] of random m x m matrices Z and can be 
evaluated by the standard tools of the random matrix theory. 

If N < m + n, i.e. N — n < m, then, generically, the matrix P*P in (37) has rank N — n 
and, as a consequence, z = 1 is an eigenvalue of Q*Q of multiplicity m + n — N . Therefore, in 
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this case the image of U (N) under the map w is a submanifold of the boundary of the matrix 
ball Q*Q < I m . This submanifold is defined by the equations 



d k 

— det(zI m -Q*Q) 



= 0, fc = 0,l,...,m + n-JV-l. (40) 

2 = 1 



An expression for dpN,nxm in this case can be obtained by following the method of calculation 
of dpN t n Xm for N > n + m that was suggested in [13] and extended in [11]. It makes use of the 
matrix integral [14] 

/ = = Q k e izTrX detX l - k , !>fc,lniz>0andl>0 

J det l {F-zI k ) ' ~ 

over Hermitian k x k matrices F. Extension of this calculation to the case of N < n + m 
requires evaluation of the above integral for positive semi-definite matrices X of rank I < k. 
Such evaluation was given in [23]. Using this result one can calculate dpN,nxm for N < n + m, 
however, the resulting expression contains many delta-functional factors and is not very useful 
for direct applications. Fortunately, evaluation of integrals (35) is a simpler task and can be 
accomplished with the help of the following standard calculation from multivariate analysis. 

Consider C lxk , the space of complex I x k matrices. If I > k then any matrix P G C ixfe of 
rank k can be uniquely written as P = HT where T = (T^) G <C kxk is upper-triangular with 
positive diagonal elements and H E C lxk is such that H*H = Ik. Correspondingly, the cartesian 
volume element (dP) in C lxk transforms as follows, see e.g. [29, 28], 

k 

(dP) = const x det(TT*y- fc ]J T^ k ~ i)+1 (dT){H*dH) (41) 

i=i 

where 

(dT) = j | Tu Q dReTijdlmTij 

l<i<k l<i<j<k 

and (H*dH) is the invariant volume element 1 on the Stiefel manifold Vk(C l ) = U(l)/U(l — k) of 
complex I x k matrices with orthonormal columns. It follows from (41) that if / is a function 
on C kxk then 

/ f(P*P)(dP) = const, f f(Z*Z)det(Z*Z) l - k (dZ), (42) 
Jcixk Jc kxk 

where we have used (41) twice, at first making the substitution P = HT and then making the 
reverse substitution T = V~ l Z with Z G £ kxk an d V G U(k). Similarly, if g is a function on 
C lxl then 

/ g{PP*){dP) = const. / / g{HZZ*H*)det{ZZ*) l ' k {dZ){H*dH). (43) 
Jc lxk Jv k (C) Jc kxk 

The integration rules (42) and (43) hold for I > k. If I < k then 

/ f{P*P){dP) = const, f [ f(HZ*ZH*)det(Z*Z) k -\dZ)(H*dH) (44) 

JC lxk JVi(C k ) Jc lxl 



1 For details of its construction, see [22, 29] 
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and 

/ g{PP*){dP) = const. / g(ZZ*) det(Z Z*) k ~ l (dZ) . (45) 

j0Xk JC lXl 

Returning to integrals (35), let us consider, alongside with uj, another map defined by par- 
tition (36), 

It maps the unitary group U (N) onto the Stiefel manifold of N x m matrices H with orthonormal 
columns, H*H = I m . Obviously, the image of the Haar measure dfj,n(U) under this map 
is invariant with respect to the right and left multiplications by unitary matrices. Since the 
Stiefel manifold is a coset space of the unitary group, such invariant measure is unique up to a 
multiplicative constant. Therefore, 

/ 9 {t(U)) dfi H {U) = [ g(H) dii H {U) = const. / g(H)(H*dH). (46) 

JU(N) JU(N) JV m {<C N ) 

Sometimes it is convenient to write the invariant measure on V^C^) as a singular measure on 

/ g(H)(H*dH) = const, f g(H)S(H* H — I m ){dH), (47) 
JV m (C N ) Jc Nxm 

where (dH) is the cartesian volume element in C Nxm and 5{H*H — I m ) is the matrix delta 
function on the space of Hermitian matrices, 

m 

5(A) = J] S(Ajj) J] <5(Re A jk )S(Tm A jk ). 

j=l l<j<k<m 

Thinking of the truncation u : U —> Q, where U G U(N) and Q G C nxm , as a composition 
of the two successive truncations 



'< : - [ p ) ^Q- 



we have, by (46)-(47), 



/ f(Q*Q)d PN , nxm (Q) = const. / ( f 5(Q*Q + P*P - I m )(dP)) f(Q*Q)(dQ). 

JQ*Q<I m Jc nxm \JC( JV - n ) xm / 

(48) 

If N > n + m then, by making use of the integration rule (42), one can replace the integration 
over Q and P in the integral on the right-hand side by integrations over m x m matrices Z and 
F, respectively, thus reducing the integral to the following one 

f ( f 5(Z*Z + F*F - I m ) det{F*F) N ~ n ~ m (dF)] f(Z*Z) det(Z*Z) n - m (dZ) . 

Performing the integration over F one obtains 

/ f(Q*Q)d PN , nxm (Q) = const. / f(Z*Z)det(I m -Z*Z) N ~ n - m det(Z*Z) n - m (dZ), 

JQ*Q<Im JZ*Z<I m 

(49) 
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in agreement with (39). 

If N < n + m, i.e. N — n < m, then by making use of the integration rule (42) to replace 
integration over Q by integration over Z as above and the integration rule (44) to replace 
integration over P by integration over FH*, where H G T^v_„(C m ) and F £ £( N ~ n ) x ( N - n ) ; one 
arrives, after carrying out the integration over Z, at 

/ f(Q*Q)dpN,nxm(Q) = const. X 

JQ*Q<l m 

f f f(I m - HF*FH*) det(F*F) n+rn - N det(I N „ n - F*F) n ~ m (dF)(H*dH). 

JV N - n (C™) JF*F<I N _ n 

If the function / is invariant with respect to the conjugation by unitary matrices, i.e. f(UAU*) = 
f(A) for unitary U, then 



f(I m - HF*FH*) = f I I m - ( F * F ° ) ) for any H such that H*H = /, 





and the integral on the right-hand side simplifies, 
f{Q*Q)dp N)nxm (Q) = const, x 



Ic 



Q*Q<In 



! f( Z 7 r ° )det(I N . n -Z*Zr +m - N det(Z*Z) n - m (dZ) (50) 

JZ*Z<I N _ n V U 1 m+n-N J 

If the function f(A) is invariant under the conjugation by unitary matrices then it is effec- 
tively a function of the eigenvalues ai, . . . ,a m of A = (^4jfc)™ fc=1 . There is one class of such 
functions for which the integral in (35) can be easily computed, see (17). These are the Schur 
functions s\ (5). The Schur functions are the characters of irreducible polynomial representa- 
tions of the general linear group. Such representations remain irreducible when restricted to 
the unitary subgroup of the general linear group. The orthogonality of matrix elements of irre- 
ducible representations as functions on the unitary group then implies the following integration 
formulae, see e.g. [27], p. 445, 

/ s x (AUBU*)dp H (U) = S ^ A >^ (51) 

JU(m) Sx(l m ) 

and 

/ s x {AU)s^BU)dp H {U) = V SA( ^? - (52) 

JU{m) Sx{l m ) 

Our calculation of the integral in (35) for / = s\ rests on the integration formula (51) and the 
invariance of the measure dpN,nxm with respect to the left and right multiplications by unitary 
matrices. 

Let X be an n x n matrix and Y be an m x m matrix. Consider the integral 

/ s x (AUBU*) dp H (U) = SA( ^ A f } (53) 

JU(N) S\{1 N ) 
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where A and B are the block diagonal matrices A = diag(X, 0), and B = diag(Y, 0). It is 
apparent from (5) that 

< n n\ J sx(a 1 ,a 2 ,...,a m ) if l(X) < m . . 

s A (ai,o 2 ,...,a m ,0,...,0) = j Q if Z(A) > m ^ ' 

where /(A) is the length of A (the number of non-zero parts A-,) 2 . Therefore the integral in (53) 
vanishes for /(A) > m (recall that m < n) and is equal to s\(X)s\(Y) / s\(In) if i(A) < m. 

On the other hand, the non-zero eigenvalues of the N x N matrix AUBU* coincide with 
those of the mxm matrix XQYQ* where Q is the principal nxm block of U, see (36). Therefore 

/ s x (AUBU*) dn H (U) = f s x (XQYQ*) dp N , nxm (Q). (55) 

JU(N) JQ*Q<I m 

The measure dpN,nxm(Q) is invariant with respect to the right multiplication by unitary matri- 
ces. Hence 

I s x {XQYQ*) d PN , nxm (Q) = [ [ s x {Q*XQVYV*) dp N , nxm (Q)dp H (V), 

JQ*Q<I m JU(m) JQ*Q<I m 

where we have used the invariance of s\(MiM2 ■ ■ ■) under the cyclic permutations of the matrices 
Mj. Reverting the order of integrations and applying the integration rule (51), one obtains 

/ s x (XQYQ*) dp N , nxm (Q) = f s x (XQQ*) dp N , nxm (Q). 

JQ*Q<I m S\{l m ) jQ*Q<I m 

The measure dpN,nxm(Q) is also invariant with respect to the left multiplications by unitary 
matrices. Repeating the above procedure, one decouples X and QQ* thus obtaining 

/ sx(XQYQ*) d PN , nxm (Q) = S j?\ Sxi ^ T \ I s x (Q*Q) dp N , nxm (Q). 

JQ*Q<I m Sx{lm)S\{l n ) jQ*Q<I m 

On comparing this to (53) and (55), one concludes that 

/ s x (XQYQ*) dp N , nxm (Q) = s ^ X >^p (56) 

JQ*Q<I m s x {1n) 

of which (17) is a particular case of X = I n , Y = I m . It is worth mentioning that the quotient 
s\(I n )s\(I m )/s\(lN) can be easily evaluated in terms of A/s by recalling Weyl's dimension 
formula 



r 

S\{In) = < 



m 



rr (\t-i-x xTT ( " + \,/ )! „ (57) 



l<i<j<m 



which holds for any integer n > m > l(X). 

By repeating the argument which was used to evaluate the integral in (56), one can extend 
the integration formulae (51) and (52) from integrals over unitary group to integrals over complex 
matrices provided that the integration measure is invariant with respect to the left and right 



2 By convention, one does not distinguish between partitions which differ merely by the number of zero parts, 
i.e. (Ai, . . . , A m ) = (Ai, . . . , A m , 0, . . . , 0) 
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multiplication by unitary matrices. For example, if L and M are mx n matrices and dp(Q) is a 
measure on C nxm invariant with respect to the right and left multiplication by unitary matrices 
and such that the integral below converges absolutely then 



/ 



s \{LQ)s^(MQ) dp{Q) = S 4^-J7 \ I s x (Q*Q)dp(Q). (58) 



In particular, if dp is the projection of the Haar measure on the matrix ball Q*Q < I m , i.e., 
dp = dpN,nxm, then (58) and (56) imply the orthogonality relation (19). 

We would like to make two remarks at this point. One is that equations (49) and (50) 
effectively give the joint probability distribution of the singular values of the truncations of 
random unitary matrices. This distribution has an interesting symmetry: Consider two square 
truncations, Q\ of size m x m and Q2 of size (N — m) x (N — m). Assuming that 2m < N, 
x = 1 is the singular value of Q2 of multiplicity N — 2m and the remaining m singular values of 
Q2 have the same distribution as the m singular values of Q\. 

The other is that some of our calculations can be repeated, almost verbatim, for truncations 
of random orthogonal matrices. In particular, let Q be the top left block of size nxm of random 
orthogonal matrix O of size N x N, m < n < N. Denote by dp (N),nxm t ne image of the Haar 
measure under the map O i-> Q. Then, by repeating the steps of the above derivation of (49) 
one obtains the integration formula 



Sc 



f(Q t Q)dp 0{N):nxm (Q) = const. 

/ f(X*X) det(I m - A^X)^-™-"^ 1 ) det(X*X)l("- m )(dX), N > n + m, 

where the integration on the right is over real m x m matrices X. This integration formula is 
not new and was previously obtained in [11] by a different method. If N < n + m and f(A) 
is a function on ^ mxm which is invariant under the conjugation by orthogonal matrices then a 
similar formula holds 



Ic 



f(Q*Q)dp ( N ),nxm(Q) = COnst.X 

Q l Q<I m 

I f( X n X t ° )det(I N - n -X*X)^ n+m ' N+1 ^det(X t X)^ n - m ^ (dX). 

JXtX<I N _ n \ U Im+n-N J 

Formula (51) also has its analogue for orthogonal matrices 

C x {XOYO l ) dp H (0) 



I 

JO(m) 



f ^_Cx(X)C x {Y) 



C x (I m ) 



where C\ are the so-called zonal polynomials. For the definition of zonal polynomials and their 
properties see [29]. As a consequence, formula (56) also has its analogue for real matrices: 
one just replaces Schur functions in (56) by zonal polynomials C\. However, the orthogonality 
relations (52) and (58) do not seem to have analogues for real matrices. 

Now we are in a position to derive the bosonic CFT formula in the range < m < N. Our 
approach is based on the Schur function expansion for the exponential e Tr M combined with the 
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orthogonality relation (58). This yields the bCFT in the form (9). The integration formulas 
(49) and (50) then lead to the specialization (2) in the interval 2m < N and (11) in the interval 
N <2m< 2N. 

Recall that the matrices X and Y are N x m. The singular value decomposition for XY* 
reads 

XY* = VDW*, D = ( p ° V (59) 

where V and W are unitary matrices of size N x N and D is a diagonal matrix of size mxm. 
The entries of D are exactly the square roots of the eigenvalues of the matrix X*XY*Y, so that 

s x (D 2 ) = s x {X*XY*Y) (60) 

In view of (59) the left-hand side of (9) takes the form of an integral over truncations of unitary 
matrices: 

/ eW Y * u + u * YX *U» H (U) = J e^ u+u ^dp H (U) = J e Tr ^ Q+ ^W xm (Q), 
U(N) U(N) Q*Q<I m 

(61) 

where Q is the top left mxm block of U and dpiy,mxm(Q) is the image of the Haar measure 
under the map U i-> Q. 

The exponential function e TrM is a symmetric function of the eigenvalues of M and as such 
can be expanded in Schur functions, see e.g. [4]: 



e TrM 



J> A a A (M). (62) 



The sum above is over all partitions A of length l(X) < N where N is the matrix dimension 
of M. However, if the matrix M has only m < N non-zero eigenvalues then, in view of (54), 
the sum is effectively over partitions A of length l(X) < m. The coefficients c x in (62) can be 
computed in terms of A: 

1 ,Y = rWA-«-A, + j), m 
\(*i-3+'V) IJ=I n;.i((+^-j)! 

where I is the length of partition A and, by convention, = for negative k. Note that c x are 
independent of the matrix dimension N. In Appendix we give an analytic method of computing 
the coefficients of Schur function expansion for multiplicative functionals of eigenvalues of M, 
to complement the algebraic method of [4]. 

Let us now expand e TrZ)< ^ and e Tr ®* D in the integral on the left in (61) in the Schur functions 
and apply the orthogonality relation (58) (recall that dpN >mxm (Q) is invariant with respect to 
the right and left multiplications by unitary matrices). This yields 

f e^ D ^ D U PN , mxm (Q) =Y,c\ 8x{X l™* Y) [ s x (Q*Q)dp N , mxm (Q), 

JQ*Q<I m x S^m) JQ*Q<I m 

(64) 
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where we have used (60). Applying (58) again, now in the opposite direction, and folding the 
Schur function expansions for the exponential functions, one obtains 

E c i Sxi ™* Y) f s,(Q*Q)d PN , mxm (Q) = [ e^™^ Y ^dp N , mxm (Q), 
A Sx^m) JQ*Q<I m JQ*Q<I m 

(65) 

and hence the bosonic CFT for m < N in the form of equation (9). 

The right-hand side of (9) takes different forms depending on m. If N > 2m then dpN^ mX m{Q) 
const. det(7 m — Q*Q) N ~ 2m , see (38), and we are back to formula (8). In the range N < 2m < 2N 
the measure dpN,mxm{Q) is n °t very handy in its explicit form. However, by (65), the integral 
on the right-hand side in (9), is effectively written in terms of the integrals 

/ sx{Q*Q)d P N ,mxm (Q), (66) 

JQ*Q<Im 

which makes it possible to evaluate it explicitly by making use of (66). 

Recall, that in the range N < 2m < 2N the measure dpN^ mxm is supported by the set T of 
mxm matrices Q such that the rank of I m — Q*Q is N — m and Q*Q < I m . This set can be 
parameterized by matrices Quzv*i see (10), and the integral in (66) is effectively over matrices 
Z, see (50). Such parametrization allows one to rewrite the right-hand side of (9) in a more 
explicit form. Indeed, it is apparent that 



Z* Z 

s a ( n j ) = sx(Quzv*Quzv*)- 

U ilm-N ' 



Hence, by (50), 



/ s x (Q*Q) dp N , m xm(Q) = / s x (Q*uzv*Quzv*) dpff yN _ m (Z)dpH(U)dp H (V), 

JQ*Q<Im JU(m) 2 JZ*Z<I N _ m 

(67) 

where dp^ N _ m {Z) is the measure defined in (12). Substituting the obtained expression for 
integral (66) into (64), one obtains the variant of the CFT in the range N < 2m < 2N as 
presented in (11). 

We would like to finish this section with two observations. Firstly, instead of reducing the 
group integral on the left-hand side in (9) to an integral over truncations of unitary matrices, one 
can directly expand the exponentials in the group integral and then apply integration formulas 
(52) and (17), the latter in the right-to-left direction. On this way one easily obtains the bCFT 
in the form of equation (15). 

Secondly, our derivation of the bCFT does not use the explicit expression (63) for the coef- 
ficients cx of the Schur function expansion for the exponential function. All that is needed of 
ca's is the property 

C(Ai,...,A m ,0,...,0) = c (Ai,...,A m )- (68) 

Therefore our calculation will go through for any convergent series g(A) = Y^x c a<Sa(^4) provided 
that the coefficients cx satisfy (68) yielding the following generalization of the bCFT formula 
(9) 

\g{(XY*U)\ 2 dp H {U) = [ g(X*XQ)g(Y*YQ) dp N , mxm (Q). (69) 

JQQ*<I m 



j 

JU(N) 
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4 CFT, reproducing kernels and Selberg integrals 

In this section we investigate relations between the CFTs and several interesting matrix integrals. 

It is apparent from the calculations in the previous section that the bosonic CFT (9) implies 
the identity (16) and vice versa. 

Consider now the matrix integral 



d ^(U) = V- Sjjm) (B * A) m) 

u(N) det(I N -AU) m det(I N -U*B*) m ^ s x (I N ) y ' 

where A* A < Ijy and B*B < 1^. One can see that the integral on the left-hand side coincides 
with the series on the right-hand side by recalling the Cauchy identity (6). It is a natural 
generalization of the well-known expansion of the inverse determinant 1/ det(J — M) in terms 
of the complete symmetric functions h r of the eigenvalues (z\, . . . , z m ) of M, 



1 



m | x 

II 1 _ . = XI h r(zi,---,Z m ). 



det(J-M) . , _ ; o 

Expanding each of the two determinants in (70) with the help of (6) and then applying the 
orthogonality relation (52) one obtains the series on the right-hand side in (70). By making use 
of (16) and (6), one can fold this series back to a matrix integral, now over matrices Q: 

<Ihh{U) f dp N)inxm (Q 



[ dp N , mxm (Q) < N 

Jo*Q<i m det(/ mJV - Q*Q ® B*A) ' - " v ' 



! U(N) det(I N - AU) m det(/jv - U*B*) m J Q * Q <i m det(/ mJV - Q*Q ® B*A) ' 

Thus, the bosonic CFT (9) implies (71) and vice versa. 

In the range 2m < N the measure dpN, m xm(Q) coincides with dpff m (Q) of (1), and (71) 

J U{N) det(I N - AU) m det(I N - U*B*) m J Q *Q< Im det(I - Q*Q B*A) ' 

This identity was obtained in our earlier work [16]. 

In the range N < 2m < 2N one can replace the integration over matrices Q by integration 
over matrices Quzv*, see (10) and (67), and (71) takes this form 

dp H {U) f dn%N_ m {Z) 



f dp H (u) r 

J det{I N - AU) m det(I N -U*B*) m J 



det(J - Z*Z ® B*A) det(I N - B*A) 2m - N ' 

U(N) Z*Z<I N ^ m 

The integration on the right-hand side is over complex (N — m) x (N — m) matrices Z and the 
integration measure is given in (12). This identity is new. 

If m = N then dpN,mxm is the Haar measure on U(N), and (71) reads 

dp H (U) 1 



l u{N) det(7jv - AU) N det(I N - U*B*) N det(I N - B*A) N ' 
This identity is almost apparent in view of the orthogonality relation (52). 
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In the previous section we verified identity (17), of which (16) is a special case (m = n). 
Identity (17) implies 

dfi H (U) f dp N)nxm (Q) 

m,n < N. 



f dfj, H (u) = r 

J U(N) det(I N - AU) m det(I N - U*B*Y Jo 



j u(N) det(/jv - AU) m det(I N - U*B*) n J Q *Q<i m det(I mN - Q*Q ® B*A) ' 

(72) 

in the same way as (16) implies (71). In the range m + n < N one can use the explicit expression 
(38) for dpN t nxm(Q)- Replacing the integration over the n x m matrices Q by integration over 
k x k matrices Z, k = min(m,n), as in (39), one obtains the identity which was claimed in (20). 
Consider now the matrix integral 

/ det(I N + AU) m det(I N + U*B*) m dp H (U) = V s y( B * A ) ( 73 ) 

Ju(N) *y s\'Vn) 

where A' is the partition conjugate to A. This integral is dual to the one in (70). The equality 
in (73) is a straightforward consequence of the orthogonality relations (52) and the dual Cauchy 
identity (7) which is a generalization of the well known expansion of the determinant det(/ + M) 
in elementary symmetric functions e r of the eigenvalues zi, . . . , z m of M, 

m m 

det(/ + M) = Y[(l + Zj ) = J2 e r(zi,...,z m ). 

j=l r=0 

In [16] we proved identity (21) which is dual to (16). It is apparent that (21) implies the matrix 
integral claimed in (22) and vice versa. We shall now give an independent derivation of (22) 
from the fermionic version (3) of the CFT and thus showing that the fCFT implies (21). 

Recall that detM can be written as a Gaussian integral over anticommuting variables tpj, 
(fj with Berezin's integration rules f (1, (fj)dipj = (0, 1), f (1, ip*)d<p* = (0, 1), 

detM = j e^^^^Yldifjdip* = J e <"* M <" (dtp). 

j 

By doubling the dimension, 

det(/iv + AU) det(Jjv + U*B*) = - det 
and 

det(/jv + AU) m det(I N + U*B*) m = ±JJ e T ^*( u+B ^ + **( u * +A W(d<5>)(d*). (74) 

Here <3? and ^ are N x m and <I>* and * are m x N matrices with anticommuting entries and 
(d&) (correspondingly (d^f)) is the product of the "differentials" of the matrix entries of $ and 

(correspondingly, d 1 ^ and d^*). The sign in front of the integral in (74) depends on the 
particular ordering of terms in these products. It is not essential for our calculation (one can 
always verify the right sign at the end of calculation) and will be omitted. 

On substituting (74) in the integral on the left-hand side in (73) and applying the fCFT (3) 
one reduces this integral to the following one 

(((I e^' B '* + ** A *+^-Q'*'^dp F Nm (Q)\ (d$)(d<Z) (75) 
J J \Jc mxm ' / 



U* + A 
U + B* 
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The quadratic form in the exponential, 

Tr (tf + <$>*A^ + *tf Q - Q*$*$) = 

m m 

E WiVvj + V^) + E (Q^V* - . 

jr' = l «J=1 

where ^ - and • are the columns of $ and and ip*j and V*j are the rows of <E>* and ^>*, is 
the one defined by the matrix 



M 



-Q* ®In I m ®A 
I m ®B* Q®I N 



Therefore, / / . . . in (75) yields det M = det(Q*Q I N + I m ® B*A) and one arrives 

at the identity 

f det(I N +AU) m det(I N + U*B*) m dfi H (U) = [ det(Q*Q®I N +I m ® B*A)d^ F N m (Q). (76) 

JU(N) J C mxm 

On making the substitution 

Q = Z-\ {dQ) = det(Z*Z) 2m (dZ), 

in the integral on the right-hand side in (76), one obtains the matrix integral (22). 

Identity (16) and its dual version (21) are rather useful in the context of Schur function 
expansions. For example, the matrix integrals (25) and (26) are straightforward corollaries 
of these identities. Consider, for example, the matrix integral in (25). Expanding each of the 
determinants on the left-hand side with the help of (6), one arrive by the way of the orthogonality 
relation (58) at 

/ T777 y d f^NA Q } (T TT^W = E 4tt4 s x {Z x Z^)f Sx (Q*Q)d PN , mxm (Q)- 

jQ*Q<i m det{I m - Z X Q *) N det(I m -QZ*) N ^ s\{l m ) J Q *Q<i m 

In view of (16), the series on the right-hand side folds to 

? SX{IN) Sx{ZlZl) = det(I m -Z lZ *r- 

A 

Hence (16) implies (25) and vice versa. Similarly, (21) implies (26) and vice versa. 
If N > 2m then dp N . mxm — df^N m and (25) reads 



■Jc 



i (77) 



*\N ' 



Q , Q<Im det(7 m - Z 1 Q*)^ det(J m - QZ\Y det(/ m - Z X Z\) 

The matrix integrals (77) and (26) are variants of integrals obtained by Berezin in [5, 6]. We 
would like to elaborate more on this link and quote some of Berezin's results. Consider complex 
rectangular matrices Z and define 

fl B = {Z £ C nxm : Z*Z<I m } , n F = C nxm , (78) 



19 



and, (cf. (1) and (4)) 

d^N,nxmi z ) = const - det(J - Z *Z)~ n ~ m+N (dZ), Zen B , N>n + m, 
d^N,nxm( z ) = const. det(J + z*ZY n ' m ' N (dZ), Zett F , N > 0, 

where (dZ) is the cartesian volume element in C nxm . The multiplicative constants are fixed by 
the normalisation j n dfiN,nxm(Z) = 1. 

Let Ab be the Hilbert space of analytic (in the nm variables Zij) functions on Qb with the 
scalar product 

(/,<?)= / j{ZUZ) d^ >nxm (Z) (79) 

JQ, B 

and Af be the Hilbert space of analytic (in the nm variables Z^) functions on Hp with the 
scalar product 

(/,<?) = / f{ZMZj d^ nxm (Z). (80) 

For (/, /) in (80) to be finite, f(Z) cannot grow at infinity faster than a certain power of ||Z||. 
Hence Af consists of polynomials in Z^ and is a finite-dimensional subspace of L 2 (Qf, d^ nxm ). 
Let us choose a basis fj(Z) in Af- Given a vector / in I? (ftp , du^ nxm ) , its orthogonal projec- 
tion on Af is given by 

(K F f)(Z) = Y,(fJj)fj(Z)= I K F (Z,Q)f(Q)d^ nxm (Q), 

where Kf(Z, Q) = fj(Z)fj(Q). As Kp{Z, Q) is the kernel of the operator of orthogonal pro- 
jection onto Af, it is independent of the choice of basis there. Similarly, if fj is a basis in Ab then 
Kb(Z, Q) = J2j fj{Z)fj{Q) is the kernel of the orthogonal projection from L 2 (^B,dfiff nxm ) 
onto Ab- The Hilbert space Ab is infinite dimensional and, hence, there is a question of con- 
vergence of the series for Kb(Z,Q). It can be shown that this series converges absolutely and 
uniformly on compact sets, see [5]. The kernels Ki(Z,Q) define coherent states in Ai, i = B,F. 
Indeed, for each Q £ Oj, Ki(Z, Q) as a function in Z, /q(Z) = K{(Z, Q), belongs to Ai and for 
any / 

= £(/>£)/;■(<?) = /(<?)■ ( 81 ) 

j 

Hence, one has the resolution of identity 

(/,<?)=/ (f,fQ)(fQ,9)dt4 f , nxm (Q), i = B,F. 

The vectors fq are not orthogonal. Putting f% for / in (81), one obtains (/^,/q) = fziQ)-, or > 
equivalently 

/ K l (Q,R)K l (R,Z)d f i i N , nxm (R) = K l (Q,Z) i = B,F, (82) 
Jfii 

which means that the kernels Ki(Z,Q) are reproducing. The coherent states which are defined 
by these kernels played an important role in Berezin's construction of quantization in symmetric 
spaces. For some parameter values these kernels can be found in the explicit form [6] 

K B (Z, Q) = det(/ m - Q*Z)~ N for real N>n + m, 
K F (Z,Q) = det(I m -Q*Z) N for TV = 0,1,2,..., 
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and identity (82) expressing the reproducing property of these kernels takes the form of the 
matrix integrals 



/ 

Jn. 



d »N,nxm(Q) 



(83) 



/n B det(/ m - Q*Z 1 )X det(I m - Z*Q)" det(/ m - Z\Z{f 

and 

det(/ m + Q*Zi) N det(/ m + Z* 2 Q) N d^ nxm (Q) = det(/ m + Z^) N , (84) 



/ 

Jsi, 



of which (25) and (26) are particular cases. Correspondingly, the matrix integrals (83) and (84) 
allows one to extend identities (16) and (21) to rectangular matrices as claimed in (17) and (23). 
Indeed, assuming that N is a non-negative integer and expanding the integrands in (83) and 
(84) in Schur functions by the means of the Cauchy identities (6)— (7), and then applying the 
orthogonality relation (cf. (58)) 

I S X (AQ)S^Q*B*) rf/4,nxm(Q) = V TT^Tt \ l S *(Q*Q) d »N,nxm(Q), * = B,F, 
JCli S\{l n )S\{l m ) Jq. 

(85) 

one obtains Schur function expansions for the determinantal powers on the right-hand side in 
(83) and (84). On comparing the coefficients in these expansions with those in (6)-(7) one 
arrives at (17) and (23). 

Identities (17) and (23) can also be derived in an elementary way and independently of (83) 
and (84) and then used to prove (83) and (84). We shall demonstrate this at the end of this 
section. 

With identity (23) in hand, we can revisit the matrix integral (73) and consider unequal 
powers of determinants 



/ 



det(I N + AU) m det(lAT + U*B*) n d m {U) = V s ^ Im > x [ 1 ^ Sy (B*A). (86) 

, s\'Hn) 



'U(N) x 

By making use of (23), one can fold the sum on the right, 

£ Sx[Im) ( S T X \ In) sy(B*A) = [ det(J + Z*Z B* A) d^ nxm (Z). 

Replacing the integration over the n x m matrices Z by integration over k x k matrices Z, 
k = min(m,n), as in (39), one obtains the identity claimed in (24). 

Now we turn to identities (17) and (23). The integration on the left-hand side in (17) and 
(23) goes effectively over the eigenvalues of Q*Q. Consider N > 2m assuming without loss of 
generality that m < n. Recalling the singular value decomposition for Q, Q = H\fXV* , where 
V is m x m unitary, i.e. V G U(m), H is n x m unitary, i.e. H £ V m (C n ) (see Section 2), and X 
is diagonal m x m matrix of the eigenvalues x±, . . . , x m of Q*Q, one can make the substitution 
Q = H\fXV* in the integrals in (17) and (23). The corresponding Jacobian is well known, see 
e.g. [29, 28], 

m m 

(dQ) = const. T] { Xi - Xj ) 2 \\ xf~ m) Yl dXj(H*dH)(V*dV), 

l<i<k<m i=l j=l 
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where (H*dH) and (V*dV) are the invariant volume elements in, respectively V m (C n ) and U (m). 
This substitution reduces the matrix integral in (17) to (27) and the one in (23) to (28). The 
normalization constants are given by 



and 



<m = /■■■/ n *r™(i-* i ) jv - m - B u ■ -'-.h 2 n ( 87 ) 

® ® 3=1 ^<i<j<m j=l 

Y\ (1 + i)! (n-m + j)\ F(N - n - m + j + 1) 
f} Q T(N -m + j + 1) 



poo poo m x n-m 

<m = /•••/ n (1+x \N +m+n n ^ Hi- m 

J ° J ° j=l [ - l+X l> l<i<j< m , I 

(l+j)\(n-m + j)\T(N + j + l) 



3=1 ' l<i<j<m J=l 

m— 1 

Q r(iV + n + j + 1) 



The integrals in (87) and (88) are particular cases of the Selberg integral 

„\ m m m—1 



r( P + j7)r((7 + i 7 )r(i + (i+j)7) 



io io/ =1 J M J r(p + g +(m + J -l) 7 )r(l+ 7 ) 

poo poo m tP~ l m 

which is a multivariate generalization of the Euler beta integral 

/V' ( i-*)«-'<fa= r -EMTM Eh,). 

io 7o (i + r( P + «) ww 

Here T and B are the Gamma and Beta functions, respectively. 

The rest of this section is devoted to explicit evaluation of the integral in (27) and the one 
in (28). To the best of our knowledge, this evaluation is new. 

Let 

S\(p,q;m)= ■■■ s x (xi, . . . ,x m ) a^ _1 (l - Xjf' 1 ]J (x» - xj) 2 dx j 

® ® j=l l<i<j<m j=l 

and 

noo roo m x ^r l m 

S^(p,q;m)= /■■■/ s x ( Xl , . . . , x m ) J] ' +2(m _ 1} J] ~ x i) 2 II ^ 

In view of (5), 

s A (xi,...,x m ) | [ (xj - x 3 : f = det (z™ +A '~ J ) _ det (x™~A _. 
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By making use of the Gram identity 

■ ■ ■ J det {Fi( Xj ))™ =1 det (Gi(xj))™ =1 J[ dx t = ml det F i (x)G j (x)dxj , (89) 
we have 

Sx(P, m) = ml det ( [ x m+ P +f 3 -i-i^ _ x y-^dx] = ml det (B(m +p + fj-i, q))f j=1 , 
\Jo J i,j=l 

where we have introduced fj = m + Xj — j, and 

c m+p+f J -i-l dx 



S x {p, q; m) = ml det 



j x m- 



x \p+q+2{m-l) 



ml det (B(m +p + fj - i, q + m - fj - 2 + i))™ (90) 



, _ 2 -+- ?:v - 

In [16] we proved the following identity for binomial determinants 

det (B( Pj - i, q 3 + i))™ =1 = det (B( Pj - i, q 3 + 1))™ =1 . (91) 
By making use of this identity, 

S\{p,q;m) = mldet(B(m + p + fj-i,q + i-l))™ j=1 

= ml < 



Now, recalling the standard determinant 



±J i T(m + p + q + fj - 1) j lJ 



det(r(p J + m-0)^ =1 = n r fe) II fa-Pi), ( 92 ) 



one arrives at 



gfosnO-m. fi r??; 1 ?/?, 
i f = ^r(m + p + g + /j - 1) 



• II (fi-fj), J) '" • A, j. (93) 

l<i<j<m 



Similarly, applying (90) to the determinant on the right in (90), one obtains 
5f(p, 9 ;m) = m!det(£(m + p + / j -^ + m-/ j -l))™ =1 

T{q + m- f j - 1) 

If 

.r- 



m 



ml < 



and, by (92), 
S^ip, q', Tn) = ml < 



"rCp + ^rfa + m -/,-!) = m + Aj -j. (94) 
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We would like to finish this section with a calculation showing that (93) implies (83). As 
has been mentioned above, for integer N this claim follows from the Cauchy identity (6). The 
case of non-integer N can be handled with the help of the following identity where M is m x m: 

1 -Va*- fj )s\(M) 3 B_f i nN + \ j -j + l){rn-j)\ 

det(/ m - M)^ ~^ Px SA(/m)SA(M) ' & " 11 riN-j + lHm + Xj-j)]' (95) 

This identity is due to Hua [21]. In Appendix we derive it and its dual version with the help of 
integration over the unitary group, see examples 4 and 5. 

It can be verified from the Selberg integral (93) that for real N > n + m and n> m 

[ s x (Q*Q) d^ nxm (Q) = S -^. (96) 

JQ"Q<I m P\ 

By making use of the Schur function expansion in (95) and orthogonality relation (85), the 
left-hand side in (83) can be expanded as follows 

£ { ^ ? *f, m) ( [ s x (Q*Q) d^ nxm (Q)) s x (Z* Zl ). 

~ S\(ln) \jQ*Q<I m J 

Applying now (96) and then folding the series with the help of (95), one gets the right-hand side 
of (83). Hence (95) implies (83) for real N > n + m. 

A similar calculation shows that (94) implies (84) for integer A > 0. One only needs to 
recall the dual Cauchy identity (7). 



5 Deformed version of the bosonic CFT. 

In order to introduce the deformation of the bosonic CFT that we are going to derive, we would 
like to demonstrate an evaluation of the right-hand side of the bosonic CFT in terms of the 
eigenvalues of the matrix X*XY*Y. We only consider the range 2m < N, however, a similar 
calculation gives an answer in the range N < 2m < 2N. As we shall eventually replace the 
integration over the matrix ball which is a bounded domain by integration over a hyperbolic 
domain which is unbounded, we first rewrite (8) by changing to complex exponentials 

f e-^'VX+x'^dvHiQ) = const, f e -iTr(x* xq+q*y*y) de t(/-Q*Q) iV - 2m (dQ). 

JU(N) JQ*Q<I m 

(97) 

Denote Qa = X*X and Qr = Y*Y and let D is diagonal matrix of the square roots of the 
eigenvalues of the matrix product QaQr- Then one can always find a non-degenerate matrix T 
such that 

Q A = TDT* and Q R = {T^DT' 1 . (98) 

Such parametrisation is possible for any two positive definite matrices [17]. This can be seen by 
diagonalising the matrix Q][ 2 Q aQ]( 2 ■ Writing Q^QaQ 1 ^ 2 = UDU*,we haveT = Q R 1/2 UD^ 2 . 
It is apparent from (98) that T is defined up to the right multiplication by diagonal unitary 
matrices. 
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By making use of (98) and singular value decomposition for Q 

m 

Q = Udmg( qi ,...,q m )V, (dQ)^(U*dU)(V*dV) ]J (q 2 - q 2 ) 2 JJd((£), 

l<j<k<m j=l 

one rewrites the integral on the right-hand side in (97) as follows 

... U dq 3 q 3 (l-q)T-^ J] ~<&) 2 * 
J ° J ° j=l l<j<k<m 

[ dn H (U) f dfi H (V) exp[-iTr(UqVTDT* + (T*)- 1 DT^V* qU*] , (99) 

Ju(m) Ju(m) 



'U(m) JU(m) 

where q = diag((?i, . . . q m ). The integral in (99) was computed in [36] 



f dfi H (U) [ dfi H (V) 

JU(m) Ju(m) 



: Tr(UAVB+CV*DU*) K det (lo{x jyk )) j k=l 

'U(m) Ju(m) Wl<j<k<m( x j ~ x k)Wl<j<k<m(yj ~ Vk) 



where x 2 and yj are the eigenvalues of AB and BC, respectively, and Iq is the modified Bessel 
function of zero order. In our case AD = —Aq 2 and BC = TD 2 T~ 1 , and collecting everything 
together we obtain 

f e-^< Y * ux+x * u * Y ^d m (Q) = 

JU(N) 



const. 



f ■ ■ ■ f f[ d mj (l - q]) n - 2m [J (g? - q 2 k ) det ( J (2q 3 d k ))™ k=v 

J ° J ° 1=1 Kj<k<m 



{dj-t%)j jo j- x i<f<k< m 

l<3<k<m J — — 

where Jq{z) is the Bessel function of zero order. By making use of the Gram formula (89), 

det f / Jo(2^4)^ 2(m ~ i)+1 (l " q 2 ) N ' 2m dq\ 

I e-^ ux ^ u * Y UMQ) « — n ((fl ^ ^ (100) 

Ju(N) 11 {dj-a k ) 

l<j<k<m 

This formula can also be derived directly from the character expansion (65), see Lemma 5 in 
[16]. Such calculation is standard, and in fact the above two- fold integral over unitary group was 
evaluated in [36] using the character expansion method [4]. Formula (100) can also be extended 
to the range N < 2m < 2N. In this case one gets the Bessel function Jo and its derivatives in 
the determinant on the right-hand side in (100). 

Formula (100) allows one to obtain an alternative version of the bosonic CFT with the 
integration manifold in the right-hand side parametrized by the matrices Q\ and Q 2 defined in 
Eq.(31). Now, consider the integral 

F m {Q A ,Q R ) = J(dQidQ 2 )det(I m -Q 1 Q 2 ) N ~ 2m exp[-iTr(Q 1 Q R + Q 2 Q A )] 

= const. / .../ n ( P 2 -pt) 2 Up^-p") N - 2m ii^ 

J ~ 1 J ~ 1 l<j<k<m j=l j=l 

x f dii H {T)ex V {-i[TPT*Q R + {T*)- 1 PT- 1 Q A ]} (101) 
JGL m (C) 
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Recalling the Q A = T DTq and Q R = {T^DT^ 1 for some T G GL m (C), see (98), we can 
rewrite the integral in (101) in terms of P, D and T. The matrix To disappears because of the 
invariance of dnu(T). The resulting integral is known [17, 15] 



const, det (Ko(2ipjdk)) ■ , _, 

/ ^H(r)exp{-i [Tw+T-^cn- 1 ] p} = n fp2 L m — ^1 , 2) 

JGL m (C) lh<j<k<m\Pj ~ Pk) lh<j<k<m\ d j ~ d k> 

where Kq(z) is the Macdonald function. Note that because the matrices D and P are diagonal, 
the above integral is effectively going over the right cosets GL m (C)/U(l) x ... x U(l). On 
substituting this expression back in (101) and using the Gram formula (89), one obtains 

det (} K (2pd k )p^ m -^ +1 (l -p 2 ) N ' 2m dp\ 
F m (Q A ,Q R ) oc - 2 _ (102) 

l<j<k<m 

Recall the identity ([18], Equations 8.405 and 8.421) 

Ko(iu) = -\ [Y (\u\) + isgn(«) J (|«|)] 
where Yq is the Neumann function. It follows from this identity that 

f 1 K {2pd)p^ m ^ +1 (l-p 2 ) N - 2m dp = 2% f 1 J (2pd)p 2( - m ~ j ^ +1 {l - p 2 ) N - 2m dp. 
J-i Jo 

Therefore, the determinants in (100) and (102) differ only by a multiplicative constant, and we 
arrive at the variant of the bosonic CFT , Eq.(34). 

It is instructive to write the formula (34) for the simplest but yet non-trivial case of m = 1. 
In this case Q\ and Q 2 are just real numbers, 

Qi = \t\ 2 p, Q2 = ^, P€[-l,l], t = Re ie eC, 



the integration measure is 



2\N-2, dtdt . 2\N-2 - dR 



(dQ 1 dQ 2 )=p(l- P T dp =p(l -pT" 1 dp M, 



and (34) reads 



•1 poo 

p[-iTr(y*Ux + x*U*y)]d t x H (U) oc / 

'U(N) 

This should be compared to the original version of the bosonic CFT. For m = 1 it reads 



/ exp[-iTr(y*Ux + x*U*y)]dn H (U) oc f p(l-p 2 ) N - 2 dp f ^ exp Lip ( R\x\ 2 + 

JU(N) JO J-oo K L V K J 



[ exp[-iTr(y*Ux + x*U*y)]dn H (U) oc / dzdz{\ - \z\ 2 ) N ~ 2 exp \-i(z\x\ 2 + z\y\ 2 )] 

Ju(N) J\z\ 2 <l 



T(V) J\z\* 

oc / p(l-p 2 ) N - 2 dp M2p\x\\y\). 
Jo 
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Therefore the transition from the original version of the bosonic CFT to its deformed version 
amounts to replacing the Bessel function by an integral, 



Ji 



dR 
~R 



exp 



ip Ra 2 + 



R 



p,a,b > 0. 



We would like to demonstrate the usefulness of the deformed CFT (102) on the example of 
the matrix integral 



R e (BB* 



, U[N} det[e 2 I N + (I N - BU)(I N - U*B*)] m ' 
For e > this integral is well defined for any matrix B. By doubling matrix dimension, 



(103) 



det[e 2 / + (/ - UB)(I - U*B*)\ 



el i(U* - B) 
i(U-B*) el 



The quadratic form f*Mf, f* = (x*,y*), corresponding to the 2N x 2N matrix on the right- 
hand side is 

e(x*x + y*y) + ix*By + iy*B*x — ix*U*y — iy*Ux. 
By making use of the Gaussian integral 



_J_ = J_ f 

detM ir 2N J C 2N 



-f*Mf 



(df), 



we have 



R £ (BB*) oc f (dX) [ (dY)e^ T ^ xtx+Y ' Y ^ xtBY - iYtBtx h cf (X*X,Y*Y) 



where 



J cf (X*X,Y*Y) = [ dii H {U) 

JU(N) 



-iTr(Y*UX+X*U*Y) 



This integral is exactly the one appearing on the left-hand side in the bosonic CFT (2). However, 
if one mindlessly applies (2), one gets a diverging integral and the use of CFT for evaluation of 
integral (103) appears to be problematic. The deformed version (34), as we shall show below, is 
free of this problem. 

By making use of (34), 



R £ (BB*) oc f (dX) [ (dY)exp 

JirNxm JirNxm 



IQNxm. JC N> 

with the quadratic form in the exponential being 



-(x*,y*)M 



(104) 



e J2( x *j x 3 + VjVj) ~ i J2( x *J B yj + y*j B * x j) + u Yl (Q^jkxtxj + i (Q2)jky* k yj- 

j=l j=l j,k=l j,k=l 



The matrix Ai corresponding to this form is 

e(I m + iQi)® I N -U m ®B 

e(I m + iQ 2 ) ® In 



M 



-il m ® B* 



el + i 



Qi (g) I N -I m (g) B 
-I m (g) B* Q 2 ® I N 
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Since M. = e + i (Hermitian matrix), the Gaussian integral in (104) converges for e > and is 
equal to 

1 1 



detM det{[(eJ m + iQi)(eI m + iQ 2 )] ® I N + ® P*P} 
or, on substituting Qi = TPT*, Q 2 = (T^^PP" 1 , 
1 1 



detM det {^(TT*)- 1 + iP}[e(TT*) + iP] ®I N + I m ® B*B} ' 

Thus finally 



,2\AT-2m x . 



/" dfl H (U) = f rr . _ 2 y 

J U(N) det[eH N + (I N - BU)(I N - U*B*)r 7 h i,i ]m f = \ Pi PA Pj> 

rr f„2_ /" ( dT ) tt i 

JGL m (c) det(T*Tr |i det {[ e (T2-)-i + ,P][ £ (TT*) + zP] + 6]} 



l<j'<fc<m 

where P = diag(pi, . . . ,p m ) and 6| are the eigenvalues of B*B. 
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A Appendix 

In this appendix we would like to demonstrate a simple calculation, reminiscent of the one given 
in Section 3 in [19], to determine coefficients in the Schur function expansions for the class of 
symmetric functions g(z±, . . . , z m ) = \Yj=\ h(z m ), where h(z) is analytic in a neighbourhood of 

N = i, 

gfa,. ..,z m ) = ^2c\s\(zi, . . .,z m ). (105) 

A 

Thinking of the Zj's as of eigenvalues of unitary matrix U, we can rewrite (105) as g(U) = 
^2\C\S\(U). Then, because of the orthogonality of the Schur functions on U(m), see (52), the 
coefficients c\ are just the "Fourier" -coefficients of the function g(U): 

c x = ! g(U)7JU) dfi H (U). (106) 

JU(rn) 

The unitary matrix U can be brought to diagonal form by a unitary transformation, U = 
Ve lii> V*, where $ = diag(0i, . . . , 4> m ), < 4> < 2ir. Correspondingly, the volume element in 
U(m) transforms as follows (see e.g. [21] or [12]) 

m 

(U*dU) = const. Y\ \ ei<Pj ~ e ^ k \ 2 II ( V * dV )' 

l<j<k<m j=l 
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j,k=l 



and the integral on the right-hand side in (106) reduces to 

i ,-2-e ,-2-k '"i m 

CA = ^H2^ / ■ ■ ■ / II h V*) det (e^ (m+Afc - fc) ) det ( e *(-*>) JJ dfc, 
where we have also used (5). By the Gram identity (89), 

c x = det ^— y / l ( e ^) e -^(^- fc +i) cZ^J 

In other words, 

c A = det(a Afc _ fc+J )™ fc=1 , or = ^ ^ hie^e-^ dc/>. (107) 
We would like to give several examples. 

Example 1. Consider the function g(z±, . . . , z m ) = exp(£^ ■ zj). Then h(z) = e z . Expanding 
the exponential in the Taylor series, = l/kl, and we recover formula (63). 

Example 2. Consider the function g(z±, . . . , z m ) = YlJLi Y\k=i V(l — tk z j), as on the left-hand 
side of the Cauchy identity (6). Then 

n ^ 

H z ) = II i _ f = ^2 h r(h,---,t n )z r , 
k=l hZ r=0 

where the /i r 's are the complete symmetric functions. Thus a r = h r (ti, . . . , t n ) and c\ = 
det(/iA fc -fe+j). In view of the Jacobi-Trudi identity 

s x = det(h Xk -k+j) = det(e A / fe _ fc+j ), 

we conclude that c\ = s\(t\, . . . ,t n ), thus recovering the Cauchy identity (6). 

Example 3. Consider the function g(z±, . . . , z m ) = YYjLi Ylk=i0- + tk z j)i as on the left-hand 
side of the dual Cauchy identity (7). Then 

n n 

K z ) = YlCL + tk^ =^2e r (t 1 ,...,t n )z r , 

k=l r=0 

where the e r 's are the elementary symmetric functions. Thus a r = e r (ti, . . . , t n ), and c\ = 
det(e\ k -k+j)- The Jacobi-Trudi identity implies that c\ = s\>(ti, . . . , t n ), and we recover the 
the dual Cauchy identity (7). 

Example 4- Consider the function g(z\, . . . , z m ) = YYjL l V(l ~~ z jY \ a > 0. Then 

M*)=(i^ = E>(«)* r , ^(a) = ^±£. (108) 



r=0 

Zm i 



Thus a r — 7 r (a) in (107) and the coefficients in the Schur function expansion (105) of g(z±, . . . , 
are given by 

CA = det( 7 A fc -*+i(a))£ fc=1 . (109) 



29 



Note that one should assume \zj\ < 1 to ensure convergence in (108). Since (106) was obtained 
by integration over \zj\ = 1 there arises the question whether one can still use this formula. The 
answer is positive. For, in view of the homogeneity of the Schur functions 

s\(tz!, . . . ,tz m ) =t w s x {z 1 ,...,z m ), |A| = ^Aj, 

3 

one can always rescale Zj — ► tzj thus extending the domain for Zj to include the unit circle. 

The determinant in (109) can be evaluated in terms of A&. Introducing the notation pk = 
Xk — k, we have c\ = det(7 Pt . +J (a)). By making use of the identity 

7r+i(a) - 7 r (a) = 7 r+i (a - 1) 

and elementary operations on the columns of the determinant, 

det (7 Pfc +i(«)) = hp k +i(a),J Pk +2(a),... ,-f Pk+m (a)\ 

= l7p fc +i(o),7 Pfc +2(a- l),...,7 Pfc+m (a-l)| 

= hp k +i( a ),J Pk +2(a- l),...,7 Pfc+m (a-m + l)| = det(j Pk+j (a - j + 1)). 

Therefore, 

r(a + p fc + l) \™ _J^r(a + Pj + l 
The determinant on the right-hand side can now be easily evaluated, 



c A = det - t^J^l =<U ' PJ Met. , . 

l r(a-j + l)fe+j)!j., =1 ir(a- 3 + l) [ V(P*+j)!A- lfc =i 



/ i \ m m i ! 



Recalling that pj = Xj — j, we arrive at 



i<j<fc<m 



n ,f T r(a + A J -j + l)(m-j)! 

S AUmjll r(a _ j + 1)(m + A ._ . v . 



where the second equality follows from Weyl's dimension formula (57). Thus finally (cf. (6)) 

This identity folds for a > and \zj\ < 1. It appears in Hua's book (Theorem 1.2.5) where it is 
derived by a different method. 
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Example 5. Consider the function g(z\, . . . , z m ) = rij=i(l + z j) a i a> 0. Then 

h(z) = (1 + zf = jr 7r (a) z r , 7r (a) = r^+Ll ( 1U ) 

r=0 

and the coefficients in the Schur function expansion (105) of g(z±, . . . ,z m ) are given by c\ = 
det(7A fc -fc+j(a))™ fc=1 , with 7r (a) as defined in (111). As in the previous example, this determi- 
nant can be evaluated in terms of A&. Now for 7 r (a) we have 

7r+i(a) + 7 r (a) = 7 r+ i(a + 1). 

By making use of this identity and elementary operations on columns, 

det(7 Afc _ fe+j (a))™ fc=1 = det(7 Afc _ fc+m (a + m - j))™ k=l . 

On substituting the expression in (111) for j r (a) in the determinant on the right-hand side, 



The determinant on the right can be reduced to the Vandermonde determinant by elementary 
operations on columns, 




[f^riy n to-*), 

1 J l<j<k<m 



\T(a- Xj +j){m + Xj - j) 



n (A.-j-A.+fe) 

l<j<fc<m 



_ n \ t~t r(q + m- j + l)(m- j)! 
" s A(l m )li r(a _ Aj+j)(m + Aj _ j)! . 

Thus finally 

fio^r-EMuto a fe-fl r^-^^X-jV (I12) 

J— 1 A j— 1 

This identity, a companion to (110), holds for a > and |zj| < 1. If a is a positive integer, say 
a = N, then l/T(N — Xj + j) = for any partition A such that Ai > N + 1 so that the sum in 
(112) is finite. By a direct computation from the Jacobi-Trudi identity, 

(N + m-j)\) 



n ^^j^ t< n * - ' - ^ + *' ^ 

j=l v -' J J/ I \<]<k<m 



c\ -- 

and we recover a particular case of the dual Cauchy identity (cf. (7)) 



Q(l + z,-)^ = ^2s X '(1n)s\{z 1 ,. . .,z rn ). 

3=1 A 
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